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ABSTRACT 

We have studied solar cycle changes in the large frequency separations that can be observed in 
Birmingham Solar Oscillations Network (BiSON) data. The large frequency separation is often 
one of the first outputs from asteroseismic studies because it can help constrain stellar properties 
like mass and radius. We have used three methods for estimating the large separations: use of 
individual p-mode frequencies, computation of the autocorrelation of frequency-power spectra, 
and computation of the power spectrum of the power spectrum. The values of the large separations 
obtained by the different methods are offset from each other and have differing sensitivities to the 
realization noise. A simple model was used to predict solar cycle variations in the large separations, 
indicating that the variations are due to the well-known solar cycle changes to mode frequency. 
However, this model is only valid over a restricted frequency range. We discuss the implications of 
these results for asteroseismology. 
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Helioseismic and asteroseismic frequency-power spectra contain a 
rich array of oscillation peaks. The separation in frequency between 
acoustic (p) modes with the same harmonic degree (1) and consecu- 
five radial orders (n) is known as the large frequency separation, Av. 
Q.^ Insights into stellar structure and evolution can be obtained by de- 
termining Av because it can help constrain stell ar properties, such as 
• mass, radius, and l og g, w i th high precision (e.g. lKalUnger et alj2009l : 
^ IStelloetalJl2008L l2009al : iMiglio et al.ll2009h . Therefore, one of the 
first aims when analysing any new set of asteroseismic data of a solar- 
I type star is to determine the Av of the observable low-/ modes. 
• • ' The frequencies of high order p modes are exp ected to follow 
^ approximately the asymptotic relation llTassoulll 19801) and so Av may 
k>< be regarded as approximately constant. However, Av is dependent on 
, j both frequency and Z and the Av often quoted in literature is that ob- 
served for Z = modes. 

In this paper we describe three methods of determining the large 
separation (see Section |2}. The most intuitive way of determining Av 
is to use the individual mode frequencies (see Section l2Tb . However, 
sometimes the signal-to-noise ratio of the oscillations is insufficient to 
allow the robust extraction of individual oscillation frequencies. This 
can be the case for asteroseismic data and can be true for solar data, for 
example, when the duty cycle is low. It is therefore convenient to be 
able to obtain Av without determining the individual mode frequencies 
and so we also obtain Av using the autocorrelation of frequency-power 
spectra (Section|2j2} and using the power spectmm of the power spec- 
trum (Section|23}- 



It is well known that the Sun's p-mode frequencies vary through- 
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changes in the observed p-mode frequencies throughout the solar cy- 
cle, we can learn about solar-cycle-related processes that occur be- 
neath the Sun's surface. In Section [3] we show that solar cycle vari- 
ations in Av can be observed in Birmingham Solar Oscillations Net- 
work (BiSON) data. This could have important consequences for as- 
teroseismic studies, which use Av to determine fundamental stellar 
properties. We also show that a simple model can be used to predict 
the observed variation in Av (Section[3T). 

In Section |4] we discuss t he impact o f observational choices on 
the obtained A v. For example, (KholikOY & Hill (2008) determine the 
solar Av using GONG and MDI data and some intriguing differences 
are found between their results and the Av obtained here. A discussion 
of the implications of our results, for both astero- and helioseismology, 
is provided in Section|5] 



2 DETERMINING THE LARGE SEPARATIONS 

BiSON has now been collecting data for over 30 yr. The quality of 
the early data, however, is poor compared to more recent data because 
of poor time coverage. Here, we have analyzed the mode frequen- 
cies observed by BiSON during the last solar cycle i.e. from 1996 
April 11 to 2010 October 8. The precision with which p-mode fre- 
quencies can be determined is directly related to the length of time se- 
ries under consideration. Consequently, p-mode frequencies are often 
determined from time series whose lengths are of the order of years. 
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However, a compromise must be made here regarding the appropriate 
length of time series for study so that changes during the solar cy- 
cle can be resolved. The observations made by BiSON were divided 
into 182.5-d-long independent subsets. Over the total observation pe- 
riod the fractional duty cycles, or "fills", of the 182.5-d subsets ranged 
from 0.72 to 0.88. 

Each method of determining the large separation must be ap- 
plied over a certain range in frequency. Initially we considered the 
range 2500 < v < 3700 |,iHz. At frequencies lower than 2500 |,iHz 
there is very little change in frequency over the course of the solar 
cycle and so it is reasonable to assume that the large separation is 
also relatively constant. Above 3700 |xHz the relationship between fre- 
quency shift and activity c hanges, as the magnitude of t he solar cycle 
frequency shift decreases (■ Libbrecht & Woodardlll990l IChaplin et alj 
n~998) . Furthermore, above ~ 4100 |xHz it appears tha t modes ex- 
perience a decrease in fr e quency as activit y increases (iRonan et alj 
ll994l : IChaplin et alJl998h .lKholikov & Hilf use the range 2300 < v < 
4300 ^Hz and so in Section|4]we have repeated our analysis over this 
frequency range. We now describe in turn each of the three methods 
of determining Av. 



2.1 The fitted method 

The large separations can be obtained from the individual mode fre- 
quencies, v/_„. Estimates of the mode frequencies were extracted from 
each 182.5-d subset by fitting a modified Lorentzian model to the 
data u sing a standard likelihood maximization method jFletcher et al.l 
l2009h . The large separations are then given by 



Av,, 



V;.„+i - V,.„. 



(1) 



However, determination of the average Av using the individual mode 
frequencies is complicated by the fact that Av/_„ is correlated to Av/,„+i. 
Therefore a linear function was fitted between n and mode frequency, 
v;„, and the gradient of this linear fit giving the average large fre- 
quency separation, Av. The linear fit was weighted by the formal er- 
rors associated with the fitted frequencies. The fit was performed sep- 
arately for each / and the mean of the gradients was determined. The 
left-hand panel of Fig.[T]shows that a linear fit represents the data well. 
For the remainder of this paper we refer to this method of determin- 
ing the large separations as the "fitted method". Although modes with 
/ as high as 5 can be observed in Sun-as-a-star data these modes are 
unlikely to contribute significantly to the observed results as their am- 
plitudes are significantly smaller than modes with / < 2. Therefore, 
when using the individual mode frequencies to determine Av we have 
used information from modes with / < 2 only. We have estimated the 
average Av for alH < 2 and we have examined Av for each / individu- 
ally. 

Monte Carlo simulations were used to determine whether the un- 
certainties returned by the linear fit reflect the true uncertainties asso- 
ciated with this method of deteiTnining the large separation^400 ar- 
tificial BiSON-like 182.5-d time series were generated dchaplin et alj 
l2006h . The mode frequencies used to generate the data were shifted 
in the manner expected due to the solar cycle (see iBroomhall et alj 
12009 j. and Section [3T1 of this paper). Although the magnitude of the 
shift placed on each individual mode was dependent on mode fre- 
quency and Z, the size of the frequency shifts were based on an / = 
mode at 3000 i^iHz being shifted by 0.5 yMz. Simulations were also 
run using a smaller frequency shift but the results were found to be 
relatively insensitive to the size of the input shift. The ratio of the 
standard deviation of Av obtained from the simulations and the mean 
uncertainty on Av indicates by how much the observed uncertainties 
need to be scaled to reflect the true errors. The scaling factors were de- 
pendent on both the frequency range used to determine Av and which 



/ were used. The simulated data were given various BiSON-like win- 
dow functions whose fills represented the range of fills observed here 
i.e. 0.72-0.88. Although dependent on fifl, over the range 0.72 - 0.88 
the scaling factors were approximately constant. Therefore the scal- 
ing factor used to detemine the size of the errors was an average of 
the scaling factors obtained for the different simulated fills. When cal- 
culating Av over 2500 < v,_„ < 3700 [xHz and using < / < 2 the 
scaling factor was 1.5 i.e. the true uncertainties associated with the fit- 
ted method were 1.5 times larger than the uncertainties implied by the 
fit. 



2.2 The autocorrelation method 

In the past cross-correlations of frequency-power spectra were used 
to estimate the mean shifts in the frequencies of solar p modes. This 
approach is particularly advantageous when the data have a low S/N, 
which makes individual mode fitting difficult, because all of the fre- 
quenc y data in the range under consideration are used dReeulo et alj 
I1994|k This method is also applied in asteroseismic studies. For exam- 
ple, |Cani2anteeta D (I2010I) describe an automated method of fitting 
the autocovariance of frequency-power spectra of asteroseismic data 
in order to obtain, amongst other parameters, the large frequency sep- 
aration. 

For different frequency lags, Sv, the autocorrelation of a 
frequency-power spectrum was determined over a given frequency 
range. Let P{v) be the frequency-power spectrum, vi the lower bound 
of the frequency range under consideration (i.e. vj = 2500 |.iHz) and 
V2 the upper limit on the frequency range (i.e. vi = 3700 [.iHz). The 
autocorrelation, Ca, is then given by 



Ca 



{P[vi : vi] ■ -P[v, + (5v : V2 + Sv]) 
yjiPlvi : V2V}{P[vi +6v:v2+ 6vf}' 



(2) 



When cross-correlations were first used to determine solar- 
cycle frequency shifts a second-order polynomial was fitted to the 
log of th e cross-correlation function over ±5 piHz aro und a lag 
of uHz jRegulo et alJll994l : I Jimenez-Reves et aDll998h . However, 
iJimenez-Reves et al. I ( l200lh showed that more accurate results could 
be obtained by fitting symmetric Lorentzians to the main peak and 
peaks at the diurnal frequencies over a range of ±20 i^iHz. This ap- 
proach was based on the assumption that the mode peaks themselves 
are Lorentzian in shape but neglects mode asymmetries. Fitting the 
main peak and the diurnal sidebands works well when the duty cycle 
of the data is low. However, as the fill increases the amplitude of the 
sidebands decreases and other featu res become visible in the cross- 
correlation spectrum. IChaplin et al ]lj.2007) fitted a function based on 
7 Lorentzians: 1 for the central peak, 2 for the diurnal frequencies, and 
2 each for the Z = 0, 2 and / = 1, 3 overlapping pairs. 

A structure containing the same combination of Lorentzian peaks 
is observed in the autocorrelation at multiples of the large separation 
(see, for example, the middle panel of Fig. [T} and so by fitting the 
autocorrelation function of the frequency-power spectrum centred on 
a lag of ~ 135 i^iHz we can determine the average large separation, A v. 
We have fitted the peaks on both the positive and negative sides of zero 
lag and then taken the average of the fitted results. The fitted function, 
M(6v) took the form 



M(6y) = 



PtA [(F + <5r|,|) I2f 



^3 [Av -6y + (fc/|fc|)rfw]- + [(F + 5F|,|) Hf 



+ B, 



(3) 



where A was the amplitude of the central peak, fit was the height of 
the A^* peak relative to the central peak, F was the width of the cen- 
tral peak, 6T\^ was the increase in width of the / = 0, 2 and Z = 1, 3 
overlapping peaks, d\k\ was the frequency spacing between the central 
peak and the outer peaks and B was a background term. Here |Z:| = 
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Figure 1. Examples of the different methods of determining Av. Left-hand panel: Fitted method. The different symbols represent the fitted mode frequencies for 
different / (black crosses represent / = 0, red squares represent / = 1, and blue triangles represent / = 2). The errors associated with the fitted frequencies are smaller 
than the symbols. For example, the mean frequency eiTor for the plotted graph was 0.12|.iHz. The different lines show the Hnear fits to the observed frequencies 
(black solid line for / = 0, red dashed fine for / = 1, and blue dotted line for / = 2). Middle panel: AutocoiTelation method. The black line shows the autocon'elation 
and the red line shows the fitted model. Right-hand panel: PSPS method. The black line shows the power spectrum of the power spectrum and the red line shows the 
fitted model. 



represented the central peak, \k\ = 1 represented the / = 0, 2 overlap, 
\k\ = 2 represented the diurnal sidebands and \k\ = 3 represented the 
/ = 1, 3 overlap. The relative height, was fixed at unity for k = 
i.e. for the central peak. The increase in width, (Jr^ was fixed at |,iHz 
when 1^1 = and 2 i.e. for the central peak and the diurnal sidebands. 
Finally do was fixed at jiHz. 

Following IChaplin et al. I toot we assumed that the widths of 
the diurnal sidebands were the same as the central peak and fixed their 
frequencies at +1 1.57 [iHz. The Lorentzians for the overlapping pairs 
have widths that are wider than the central peak because of the in- 
fluence of different components of a multiplet. For example an / = 
mode, which has one component only, will overlap with the 3 visible 
(in Sun-as-a-star data) m components of an Z = 2 mode at 3 difi'erent 
frequencies. The distances of the overlapping I = 0,2 and Z = 1,3 
pairs from the central peak are known as the small separations, dij+2- 
We have used initial guess values for the Z = 0, 2 pair of cZo,2 = 8.9 ^iHz 
and for the Z = 1, 3 pair of dt 3 = 15.9 [xHz, based on the average ob- 
served spacings over the frequency range of interest. For the remainder 
of this paper we refer to this method of determining the large separa- 
tions as the "autocorrelation method". 

The underlying profile of a mode in a frequency-power spectrum 
is an asymmetric Lorentzian. Two important parameters in character- 
izing each mode Lorentzian are the height and width. We then define 
the power of the mode as the total area under the Lorentzian. The au- 
tocorrelation function is most influenced by the modes with the largest 
heights i.e. those modes around 3100 i^iHz. In efl'ect the obtained fre- 
quency separation represents a weighted average over the frequency 
range of interest. 

Although the statistics of a frequency-power spectrum is not 
Gaussian the autocorrelation function is the sum of many points and 
so the central limit theorem can be applied. Therefore, a standard least 
squares fitting can be performed. However, the points in the auto- 
correlation are highly correlated and, although this does not bias the 
fitted parameters, the formal Hessian uncertainties are an underesti- 
mate of the true uncertainties. Monte Carlo simulations, similar to 
those described in Section 12.11 were performed to estimate a scal- 
ing factor by which the Hessian uncertainties needed to be multi- 
plied to better represent the true errors of the fit. The scaling factor 
showed no systematic dependence on fill over the range examined 
here. We therefore took the scaling factor to be the mean of the sim- 
ulated scaling factors. When calculating the autocorrelation function 
over 2500 < v/_„ < 3700 the scaling factor was 3.5. 

The middle panel of Fig. [T] shows an example of the autocorre- 
lation. A prominent peak is observed and the model provides a good 
fit to the data. The fill of the plotted data is relatively high for ground- 



based observations (0.88) and so overlaps between the Z = 0, 2 and 
Z = 1, 3 modes are clearly visible, whereas the daily harmonic peaks 
are suppressed. 



2.3 The PSPS method 

An alternative way of determining 
tions that is frequently employ ed in 
(e.g. iRoxburgh & Vorontsov 



les 



Roxburgh & Vorontso' 
2009 



the large 
asteroseismic 



separa- 
stud- 



'2006; Roxburgh 200' 
Mossei- & ApDourchaux 20Q9; Jj ekker et al.. 2010 ; Mathuret 
20101 ; I Mosser & Appourchau 3 I2010I) is to calculate the power 
spectrum of the power spectrum (PS (^ PS). A significant peak 
is observed in the power spectrum of the power spectrum at half 
the large separation. This method is eq uivalent to determining the 
autocorrelation of the time series. Kholikov & Hilll (2008) used the 
autocorrelation of time series to determine the acoustic radius, T, 
which is related to the large separation, Av, by the following equation: 



Av = — . 

2T 



(4) 



The PS PS was determined, over the frequency range of in- 
terest. We have oversampled the data by a factor of 10 and this was 
done by adding zeroes to the end of the frequency range of interest. A 
peak is clearly visible at half the large separation (see the right-hand 
panel of Fig. [T](. A Gaussian was fitted to this peak to determine the 
value of half the separation. Two alternate ways of finding the cen- 
tral frequency of the peak are by interpolation and by determining the 
centroid of the peak. However, these methods were found to be more 
sensitive to the background noise and so we have instead used the re- 
sults obtained by fitting a Gaussian. For the remainder of this paper 
we refer to this method of determining Av as the "power spectrum of 
the power spectrum (PSPS) method". 

To determine whether the uncertainties from fitting a Gaussian 
to the 0.5Av peak are representative of the true uncertainties Monte 
Carlo simulations were performed that were similar to those described 
in Section im The simulations showed that the uncertainties were un- 
derestimated by a factor that was dependent on the frequency range 
under consideration. Different BiSON-like window functions were ap- 
plied to the simulated time series. As with the other two methods the 
scaling factor was independent of the fill over the range of fills ob- 
served here. When calculating the PS PS over 2500 < v,,„ < 3700 
the scaling factor was 1.8. 
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Figure 2. Large separations obtained from the 185. 5-d data sets using the fit- 
ted, autocorrelation, and PSPS methods. For the fitted method the mean Ay for 
modes with / < 2 was obtained, while the autocorrelation and PSPS methods 
contain information from all modes vis ible in Sun-as-a-star data. Scaled ver- 
sions of the frequency shifts obtained bv lFletcher et al.l feoid) and the 10.7-cm 
flux are plotted for comparison purposes. 



3 SOLAR CYCLE VARIATIONS IN Av 

Fig. |2] shows the large separations obtained from the 182. 5-d data 
sets using the fitted, autocorrelation, and PSPS methods. For the fitted 
method the mean Av for modes with / < 2 was obtained. As we have 
used Sun-as-a-star data the autocorrelation and PSPS methods con- 
tain information from all visible /. However, in practice, modes with 
/ < 2 will dominate as they are most prominent. Solar cycle variations 
are clearly visible in the fitted and autocorrelation method results. The 
PSPS method results are noisier and so it is difficult to see any solar 
cycle variations. 

There is an offset between the large separations obtained by the 
three methods (also see the top panel of Fig. [3}. The difference be- 
tween the fitted and autocorrelation methods can be explained in terms 
of the different frequency dependence of the methods combined with 
the frequency dependence of the large separations themselves. The au- 
tocorrelation method is heavily weighted towards the most prominent 
modes. However, the fitted method will be weighted towards lower- 
frequency modes since the frequencies of these modes can be obtained 
more precisely than high-frequency modes (because they have longer 
damping times and consequently narrower widths in frequency-power 
spectra). Weighting the linear fit between n and v„ ; by a combination 
of mode height and the errors associated with the frequencies reduces 
the disparity but only very marginally. The PSPS method is weighted 
by the power (and not the height) of the modes. Mode power drops off 
less rapidly with frequency than mode height and so the PSPS method 
is weighted towards higher-frequency modes than the other two meth- 
ods. Since Av increases with frequency over the range of frequencies 
considered here, the Av determined by the PSPS method is larger than 
the Av obtained by the other methods. 

The large separations have also been obtained from two sets of 
data that were simul ated for the solar Fi tting at Low Angular De- 
gree Group (FLAG) dChaplin et al.ll2006L top panel of Fig. |3}. The 
data were simulated in the time domain. In the first simulated data set 
(FLAGl) the properties of the oscillations were constant but in the 
second data set (FLAG2) the oscillation properties were modulated 
by the solar cycle. For the fitted method the mean Av for modes with 
/ < 2 was obtained while the autocorrelation and PSPS methods con- 
tain information from all visible /. The large separations determined 
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Figure 3. Top panel: Large separations observed in the FLAGl and FLAG2 
data. Bottom panel: Differences between the FLAGl and FLAG2 large sep- 
arations. For the fitted method the mean Av for modes with 1 < 2 was ob- 
tained, while the autocorrelation and PSPS methods contain information from 
all modes visible in Sun-as-a-star data. 



by the autocoiTelation and PSPS methods vary significantly more than 
the large separations determined using the fitted method. This indi- 
cates that the autocorrelation and PSPS methods are more sensitive to 
the realisation noise. 

The bottom panel of Fig.[3]shows the difference between the large 
separations obtained in the FLAG2 data set and the large separations 
obtained from the FLAGl data set. This difference uncovers the sim- 
ulated solar cycle effect, which alters the large separations by about 
0.05 [iHz. It should be remembered that this data set is only ~ 9 yr in 
length and so does not quite cover an entire solar cycle. 



3.1 Predicting variations in the large separations using solar 
cycle frequency correction techniques 

iBroomhall et al.l ( l2009ah corrected Sun-as-a-star p-mode frequencies 
for the effect of the solar cycle. Here, we have reversed the procedure 
applied to implement these corrections in order to predict the frequen- 
cies of p modes in Sun-as-a-star data at a particular leve l of activity. 
We hav e used the frequencies quoted in table 3 of Bro omhall et alj 
j2009bh as a reference set of frequencies. These frequencies were de- 
te rmined by fitting 23 y r of BiSON data using the methods described 
in iFletcher et al ] ll2009h . The raw fitted frequencies were then cor- 
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reeled, using the methods described in iBroomhall et alj ( l2009ah . to 
give the frequencies that would have been observed at the canonical 
quiet Sun activity level, and it is these frequencies th at we use here. 
The correction performed by IBroomhall et ailfcOOgal) comes in three 
parts: a linear solar cycle correction, a "devil-in-the-detail" correction, 
and a Sun-as-a-star correction. 

The linear solar cycle correction is well known and can be used 
to coiTect mode frequencies to a nominal activity level. The correc- 
tion is based on the assumption that variations in global activity in- 
dices can be used as proxies for low-/ frequency shifts. It is also as- 
sumed that the correction can be parameterised as a linear function 
of the chosen activity measure. Here we have used the 10.7-cm flux 
iTaoDine & Detracev 1990) for w hich these assumptions a re robust at 
the level of precision of the data llBroomhall et al.ll2009aL and refer- 
ences therein). The canonical quiet-Sun value of the radio flux is fixed 
from historical observations at {A{t))c = 64 x 10""^ Wm"- Hz"'. The 
size of the linear correction is also dependent on mode frequ ency and /. 
Here we have used the frequency dependence derived by Chaplin et alj 
j200lL l2004h . The / dependence is small for / < 3 but is included 
nonetheles s and occurs because of differen ces in the mode inertia 
fchristense n-Dalsgaard & Berthomieulll99lh . 

The "devil-in-the-detail" correction accounts for cross-talk be- 
tween the variations of different mode parameters over the solar cycle, 
and the distribution of activity levels over the period of observations. 
As we are using relatively short time series here, this effect will be 
negligible and so this correction was not performed. 

The Sun-as-a-star correction accounts for the fact that when us- 
ing Sun-as-a-star observations the rotationally split components of a 
mode have different visibilities. As the plane of the Sun's rotation axis 
is nearly perpendicular to the line-of-sight, only modes where I + m 
is even have significant visibility. Furthermore, estimates of the cen- 
troid frequencies are dominated by the \m\ = / components, as they 
are most prominent. The difference between the centroid and fitted 
frequencies depends on the level of activity over the period the obser- 
vations were made. When the solar activity is at a minimum the com- 
ponents are observed to be in a near-symmetrical arrangement and so 
the fitted centroid is close to the true centroid frequency. However, 
at moderate to high activity levels this is not the case as the mode 
components are not arranged symmetrically. For example, the \m\ = 2 
components of an / = 2 mode will experience a larger shift at high 
activities than the m = component. The magnitude of the observed 
asymmetry is related to the inhomogeneous distribution of the mag- 
netic activity over the solar surface and the spherical harmonic associ- 
ated with each visible m component. Therefore, the fitted frequencies 
diff^er from the true centroids b y an amount that is dependent on /. 
lAppourch aux & Chaplirj (l2007h describe how to make a Sun-as-a-star 
correction, which is determined using the so-called a coefficients that 
are found from fits for an unresolved Sun. 

Once the frequencies for a particular time series at a known ac- 
tivity level have been predicted we can use these frequencies to deter- 
mine the average large frequency separation using the fitted method. 
The predicted large separations are plotted in Fig. |4] along with the 
observed Av determined by the three methods. For the fitted method 
the mean Av for modes with / < 2 was obtained, while the autocorre- 
lation and PSPS methods contain information from all modes visible 
in Sun-as-a-star data. The agreement between the predicted large sep- 
arations and the large separations observed using the fitted method is 
very good. The good agreement between the predicted and observed 
large separations implies that the observed variation can be explained 
in terms of the observed shifts of the individual mode frequencies. 

Notice that the observed large separations show more short-term 
variability than the predicted Av in the declining phase of cycle 23 
and the minimum between cycles 23 and 24. It has long been ob- 
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Figure 4. Comparison between the predicted and observed Av. The different 
symbols represent the results of different methods (see legend) and the thick 
solid green line represents the model. For the fitted method the mean Ay for 
modes with / < 2 was obtained, while the autocorrelation and PSPS methods 
contain information from all modes visible in Sun-as-a-star data. 
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Figure 5. Comparison between the modelled and fitted method observed Av 
obtained over 2500 < v/ „ < 3700 |iHz using modes with I < 2. The symbols 
represent the observed Av and the dashed line represents the linear fits between 
the model and observed Av. The 1:1 relation has been also been plotted (the 
red solid line) to guide the eye. 



served that p-mode frequen cies respond to short-term changes in activ- 
ity (e.g lRhodes et alj2002h . Furthermore, a similar behaviour was ob- 
ser ved in the solar - cycle frequency shifts by IBroomhall et aO ( l2009bl) 
and lFletcher et al.l ( I2OI0I) . Fig. |2] shows that the quasi-biennial varia- 
tions observed in the large separations are in agreement with those ob- 
served in the frequency shifts by Broomhall et al. and Fletcher et al.1 
This implies that the shorter-term variations show some frequency de- 
pendence. However, the amplitude of the shorter-term variations in Av 
is smaller than the amplitude of the 1 1-yr variation, indicating that the 
frequency dependence of the shorter-term variation is weaker than the 
frequency de pendence of the 1 1-y r cycle. This is in agreement with 
the results of iFletcher et al.l ilOld) . 

Fig. |5| shows a direct comparison between the model and fitted 
method observed Av using modes with Z < 2. Although not shown here 
a similar analysis was performed for the autocorrelation and PSPS 
methods. Linear fits between the modelled and observed Av were per- 
formed and the gradients of the linear fits were within lo" of unity 
for all three methods. This suggests that the magnitude of the solar 
cycle effect is the same for each method of deriving Av and that all 
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1996 1998 2000 2002 2004 2006 2008 2010 
Year 

Fitted method :x — x— x/^o * "1=^ □ □ □ /=2 
Model: 1=0 l=^ 7=2 

Figure 6. Large separations observed for individual /. The results of the fitted 
method are plotted for each / < 2, as are the model predictions (see legend). 

three methods produce variations in Av that are in agreement with the 
model, even if the absolute values of Av differ (e.g. Fig.|4j. 

3.2 The / dependence of the solar cycle variations in the large 
separations 

As we are using Sun-as-a-star data the large separations cannot be 
determined for each I separately using the autocorrelation and PSPS 
methods. However, the fitted method can be used to examine the / de- 
pendence of the variation in Av. Fig.|6]shows the solar cycle variations 
in the large separations for individual /. There is an offset between the 
different large separations observed in the different /, demonstrating 
the I dependence of the large separations. Av increases with increas- 
ing / because the depth of the cavity in which the modes are trapped 
decreases and so the acoustic radius, which is inversely proportional 
to Av, also decreases with increasing /. 

The individual-/ large separations are generally well reproduced 
by the model described in Section lSTI Figs|6]and|7]show that the mag- 
nitude of the observed 11 -yr cycle variations increase slightly with /. 
For the model values the / = peak-to-trough variation is 0.06 [,iHz, 
whereas the / = 1 variation is 0.08 [.iHz and the 1 = 2 variation is 
0.09 [xHz. Solar cycle shifts in mode frequencies are dependent on / 
because the size of the perturbation is inversely related to the mode 
inertia, which decreases with increasing / (e.g. Libbrecht & Woodardl 
[l99Q) . However, this effect is negligible over the range of / considered 
here (/ < 2). As mentioned in Section ISTI there is an / dependence in 
the solar cycle frequency shifts that occurs because magnetic activity 
is inhomogenously spread across the solar surface. It is this / depen- 
dence that causes the differences with / in the amplitudes of the solar 
cycle variations in Av. 

Fig. |7] shows a direct comparison between the model and ob- 
served Av of each /. The model has been computed separately for each 
/. Although offset from each other the ranges of the axes are the same 
and so it is again possible to see that the 1 = 2 modes show more vari- 
ation throughout the solar cycle than the 1=1, and, in particular, the 
/ = modes. Linear fits between the modelled and observed Av were 
again performed and the gradients of the linear fits lie within 2cr of 
unity for all three /, with the / = 1 modes showing the best agreement. 



4 IMPACT OF OBSERVATIONAL CHOICES 

iKholikov &Hii]i2008l^ observed solar cycle variations in the acoustic 
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Figure 8. Av observed over a frequency range of 2300 — 4300 |.iHz. The differ- 
ent symbols represent the results of the different methods (see legend). Model 
predictions are also plotted. For the fitted method the mean Av for modes with 
/ < 2 was obtained, while the autocorrelation and PSPS methods contain infor- 
mation from all modes visible in Sun-as-a-star data. 



radius, which is inversely proportional to Av, in Global Oscillations 
Network Group (GONG) and Michelson Doppler Imager (MDl) d ata. 
To allow a comparison with the results o f iKhoUkov & Hilll ( l2008h we 
have repeated the above analysis using the same frequency range i.e. 
2300 < y < 4300 [^iHz. Kholikov & Hill determi ne the acoustic rad ius 
for individual /. Therefore we have compared the lKholikov 8i Hilj / = 
results with the / = Av found here. Before that we first discuss the 
impact of changing the frequency range on the /-averaged results. 

Fig. m shows the large separations found using the autocorrela- 
tion, fitted, and PSPS methods, and the model predictions. The fit- 
ted method was used to determine the average Av for modes with 
/ < 2. The change in frequency range has increased the offset between 
the results of the three methods. Changing the frequency range has 
had very little effect on the Av determined using the autocorrelation 
method because the autocorrelation method is weighted most heavily 
towards the most prominent modes, at ~ 3100 i^iHz. However, the fit- 
ted method large separations have been shifted downwards by approx- 
imately 0.2 |,iHz and the short term structure of the large separations 
has changed. The shift in the fitted method large separations arises 
because the obtained Av are weighted towards the lowest frequency 
modes in the range under consideration and modes at 2300 [iHz have a 
lower Av than modes at 2500 [xHz. The change in the short-term struc- 
ture of the fitted method large separations could be due to extra noise, 
which is introduced into the results by including the higher frequency 
modes, whose frequencies cannot be determined precisely because of 
short mode lifetimes. It could also be because low-frequency modes 
are more sensitive to realization noise, due to long mode lifetimes. An- 
other possibility is that the changes in the short-term structure could 
reflect the frequency dependence of the shorter-term signal. 

The large separations determined by the PSPS method have been 
shifted upwards by approximately 0.12 |iHz. As mentioned earlier the 
PSPS method is weighted by mode power and the power of modes in 
the range 3700 - 4300 [iHz is, on average, higher than the power of 
modes in the range 2300 - 2500 [iHz. Since modes in the frequency 
range 3700 - 4300 |,iHz have larger Av than lower-frequency modes 
the large separations determined by the PSPS method are increased. 

Also plotted in Fig. [8] are the modelled Av determined over the 
frequency range 2300 < y;_„ < 4300 [iHz and for modes with / < 2. As 
expected they maintain approximate agreement with the fitted method 
results (as opposed to the autocorrelation and PSPS results). However, 
the agreement between the modelled and fitted Av is clearly not as 
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Figure 7. Comparison between ttie modelled and observed Ay obtained using the fitted method. The results for the different / are plotted in separate panels. In each 
panel the symbols represent the observed Av and the dashed lines represent a linear fit between the modelled and observed Av. The 1 : 1 relation has been also been 
plotted in each panel (the black solid lines) to guide the eye. 
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Figure 9. Comparison between the modelled and fitted method observed Av 
obtained over 2300 < v;_„ < 4300 |xHz and / < 2. The symbols represent the 
observed Av and the dashed line represents the linear fits between the model 
and observed Av. The 1 : 1 relation has been also been plotted (the red solid line) 
to guide the eye. 
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Figure 10. Av observed for individual / over 2300 - 4300 |,iHz. Model predic- 
tions are also plotted. 



good over this frequency range as was found in Fig. |6] (over 2500 < 
v,.„ < 3700 nHz). 

Fig. [9] shows a direct comparison between the model and fitted 
method observed Ay for 2300 < y/,„ < 4300 |,iHz and / < 2. Once 
again, although not plotted here a similar analysis was performed for 
the autocorrelation and PSPS methods. The gradients of the linear fits 
all show departures from unity: the gradient of the fitted method is 
4.7cr from unity, the autocorrelation is 1 .6a- from unity, and the PSPS 
method is 4.5cr from unity. The observed Ay show more variation than 
was predicted by the model. 

Further investigation reveals that this is because the lower limit 
of the frequency range has been changed from 2500 |,iHz to 2300 |,iHz, 
which implies that the frequency dependence of the solar cycle pertur- 
bation on mode frequencies does not reflect the perturbation observed 
in BiSON data at low frequencies. 

Fig.llOlshows the /-dependence of the large separations using the 
frequency range 2300 < y„,/ < 4300 |xHz. Comparison with Fig. |6] 
shows that the Z = 1 large separations are less sensitive to the change 
in frequency range than the / = and / = 2 modes. Notice that the 
models do not represent the observed large separations as well over 
this frequency range. 

Fig. [TT] shows a direct comparison between the model and ob- 
served Ay for each /. Although the gradient of the linear fit to the / = 1 
mode results lies within 1 tr of unity the gradients of the fits to the 
/ = and 1 = 2 modes differ significantly from unity: 4.9it for / = 
and 2.2cr for 1 = 2. Again the observed Av show more variation than 
was predicted by the model. Therefore, the model does not represent 
the data for / = and / = 2 modes over this frequency range. 



iKholikov & mil j2008l) use the autocorrelation of time series 
to obtain the aco ustic radius. This method is analogous to the PSPS 
method used here. iKholikov & hIH use resolved Sun GONG and MDI 
data and so are able to isolate individual /. Here, since we are using 
Sun-as-a-star data we can only determine Av for individual I using the 
fitted method. Ther efore, in the left-hand panel of Fig. [TT] we compare 
the / = results of iKhoUkov & Hill (2008) with the / = Av deter- 
mined using the fitted method. IKholikov & Hill determined the offset 
for 36 d (GONG) and 72 d (MDI) time series. To allow a comparison 
with the results derived here we have plotted the weighted mean of the 
results of iKholikov & Hill over 5 and 3 points respectively, where the 
weights were given by the errors of the observations. The range of Av 
observed in the BiSON data is significantly larger than is observed by 
Kholikov & Hill. A linear fit between the results of Kholikov & Hil| 
and the model was performed. The gradients of both fits were signif- 
icantly shallower than unity: 17.6cr away from unity for the GONG 
results and 13.8cr for the MDI results. It should be noted t hat the un- 
certainties associated with the results of IKholikov &. Hilll were very 
small. While, to the eye, the gradients of the slopes look similar to 
unity, relative to the size of the uncertainties, the departures from unity 
are significant. 

The GONG and MDI Av found by iKhoUkov & ffiil are ofi'set 
from the BiSON results obtained here and from each other This sug- 
gests that not only do dilferent methods of determining Av result in 
different values (see e.g. Fig. [9) but th at also different data result in 
different values of A v. IKholikov & Hill were unable to explain this off- 
set. The autocorrelation method was used to determine Av for 182.5-d 
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Figure 11. Comparison between tiie modelled and observed Av obtained using the fitted method over 2300 < v/„ < 4300 |iHz. The results for the different / are 
plotted in separate panels. In each panel the symbols represent the observed Ay and the dashed lines represent a linear fit between the modelled and observed Av. The 
1 : 1 relation has been also been plotted in each panel (the black solid lines) to guide the eye. In the 1 = panel the red crosses represent the results obtained using the 
BiSON data. Also plotted are smoothed versions of the results obtained bv lKholikov & HiH llOOSj) for the GONG (pink diamonds) and MDI (orange asterisks) data. 



sets of data observed by the Global Oscillations at Low Frequencies 
(GOL F; lGabrieIetal]|l995l : Ijimenez-Reves et al.ll2003l: iGarci'a et alj 
12005') instrument onboard the SOlar and Heliospheric Observatory 
(SOHO) spacecraft. GOLF has been collecting data since 1996 and 
so we have been able to analyze data covering almost the entirety of 
solar-cycle 23, i.e., from 1996 April 1 1 to 2009 April 7. Although not 
shown here the BiSON and GOLF results were in very good agree- 
ment and the differences between the BiSON and GOLF Av were, on 
average, significantly smaller than the associated eiTor bars. It is pos- 
sible, theref ore, that the main offset between the BiSON Av and the Av 
observed bv lKholikov & Hilll occurs because of the different method- 
ologies employed. 



5 SUMMARY AND DISCUSSION 

The thr ee different methods res ulted in different values of Av. The re- 
sults of lKholikov & Hill ^20081 ) were offset from each other and from 
the BiSON data. The different values for the different methods can be 
understood in terms of the weighting with frequency of the methods 
and the frequency dependence of Av. However, this do es not explain 
the difference between the two sets of lKholikov & Hill results 

The observed offsets could be important for asteroseismic stud- 
ies, which use the determined Av to infer stellar properties such as 
radius and mass. Stello et al., C2009b.) showed that the following scal- 
ing relation holds over most of the HR diagram and results in errors 
that are probably below 1 per cent, 

Av I M/Me 

Avo ^(R/Rey' ^' 

where M is the mass of a star, R is the radius of as star and O denotes 
the values for the Sun. The farthest separated like-for-like values con- 
sidered here are the fitted method and PSPS method results for the 
frequency range 2300 < v„j < 4300 |.iHz, which are on average sep- 
arated by 0.66 [xHz. If we take the fitted method results to represent 
those of the Sun and the PSPS method results to represent those of 
another star and assume that the M/Mq = 1, equation [5] implies that 
the "other" star is 0.3 per cent larger than the Sun, which is within the 
errors implied by Stello et al. (2009b). 

The fitting method appears to be the cleanest way of determin- 
ing Av as it is less sensitive to the noise realization than the other two 
methods. However, it is not always possible to determine the frequen- 
cies of the individual modes because of poor signal-to-noise levels and 
poor fill. If this is the case, for the data examined here, the autocorre- 
lation method produced more stable results than the PSPS method. 

Solar cycle changes in Av were visible in the results of the fitted 
and autocorrelation methods. This again suggests that care is needed 



in using Av in determining properties of stars. The amplitude of the 
solar cycle effect is smaller than 0.66 [.iHz (see Fig.jJJ and, therefore, 
using the same logic as above, the solar cycle variations in the Sun 
would result in changes in stellar radii that are within error estimates. 
However, the Sun is a relatively quiet star and so it is possible that 
Av would vary more in other stars, resulting in systematic errors in 
estimated masses and radii. 

We have shown that the observed solar cycle changes in Av can 
be predicted using a simple model, which is valid over the frequency 
range 2500 < v„j < 3700 [iHz. This implies that the changes in Av 
are due to the frequency dependence of the solar cycle changes in v;.„. 
The correlation with the model is less good over the frequency range 
2300 < vij, < 4300 |xHz. This is mainly because the model under- 
estimates the size of the perturbation in Av at low frequencies, sug- 
gesting that the BiSON data do not follow the assumed frequency de- 
pendence at low frequencies. This could be because the assumed fre- 
quency dependence was derived using int ermediate- and high-degre e 
GONG data and not (low-/) BiSON data jchaplin et alj|200ll,l2004l) . 
The model also breaks down to different extents for the different /. 
This could be because the inertia scaling is not correct for the BiSON 
data at low frequencies; 

Recentlv Gar cia et alJ l l201(]h observed signatures of a stellar ac- 
tivity cycle in asteroseismic data obta ined by the Convect ion Rotation 
and Planetary Transits (CoRoT; e.g. iMichel et al.ll2008l) space mis- 
sion. With the prospect of longer asteroseismic d ata sets (~ 3.5 yr) be- 
coming available through, for example, Kepler dChaplin et al.ll20ld : 
iKoch et al. 2010, ) there will be opportunities to observe activity cycles 
in other stars. These observations will provide constraints for models 
of stellar dynamos under conditions different from those in the Sun. 
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